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Asymmetric Field Profile in Bose Glass Phase of Irradiated YBa 2 Cu307_,5: 
Loss of Interlayer Coherence around 1/3 of Matching Field 
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Magneto-optical imaging in YBa2Cu307_^ with tilted columnar defects (CD's) shows an asym- 
metric critical-state field profile. The observed hysteretic shift of the profile ridge (trough) from the 
center of the sample is explained by in-plane magnetization originated from vortex alignment along 
CD's. The extracted ratio of the in-plane to out-of-plane magnetization component has a maximum 
at 1/5 of matching field (£>$) and disappears above B$/3, suggesting a reduction of interlayer co- 
herence well bellow B<s> in the Bose glass phase. Implications are discussed in comparison with the 
vortex liquid recoupling observed in irradiated Bi2Sr2CaCu20g+y. 



An effective way to enhance the critical current density 
(J c ) in high temperature superconductors (HTSC's) is 
introduction of the correlated disorder such as columnar 
defects (CD's) 0-g and planar defects [@§. Heavy-ion 
irradiation is one of the most controllable techniques to 
introduce CD's in the sample whose diameter is of the or- 
der of coherence length in HTSC's. When the density of 
CD's, or dose-equivalent matching field B$, is increased, 
the enhancement of J c persists at higher fields (l) . How- 
ever, the field dependence of J c , which is proportional to 
the irreversible magnetization 0] , has a maximum not at 
Bq, , but at a field significantly smaller than B$> . 

Recent studies on the vortex phase diagram of highly 
anisotropic E^S^CaCr^Os+^BSCCO) with CD's have 
suggested a nearly temperature independent boundary at 
- §,§. In the vortex liquid (VL) regime above the 

irreversibility line, dramatic enhancement of the c-axis 
coherence at (1/5 ~ 1/3)5$ has been demonstrated by 
Josephson plasma resonance (JPR) studies §]. In addi- 
tion, the c-axis resistivity decreases at similar fields con- 
sistent with the enhancement of interlayer coherence [jl0| . 
Reduction of reversible magnetization p pl|p^ ] is also re- 
ported in the same field range. Theoretically a Monte 
Carlo simulation by Sugano et al. [jl3| suggests a field- 
driven transition with enhancement of vortex trapping 
rate by CD's at _B$/3. Their calculation shows that the 
interlayer coherence jumps at the transition consistent 
with JPR experiments in the VL phase of BSCCO. In 
less anisotropic YBa2Cu307_ d -(YBCO), no JPR experi- 
ments have been reported because the plasma frequency 
in YBCO is much higher than the accessible frequency 
range of JPR experiments |l . 

An interesting conclusion in Ref. Jb| is that in the 
solid phase below the irreversibility line, which is called 
as Bose glass (BG) phase [H, the sign of the interlayer 
phase coherence change becomes opposite: i.e., it de- 
creases at _B$/3. In the BG phase, the peak in J C (H) 



is actually located at (1/5 - 1/3)5$ in BSCCO |i and 
similar behavior can be found in YBCO as well [pj . In 
the collective pinning theory fl5j , the reduction of the col- 
lective pinning length related to the interlayer coherence 
causes the enhancement of the critical current p6| . Thus, 
the origin of the J C (H) peak in HTSC's may be related to 
this i?$/3 boundary. Note that in unirradiated BSCCO 
a steep magnetization increase JUT] and an abrupt reduc- 
tion of the phase coherence flg|1 occurs simultaneously 
at the second magnetization peak field. However, so far 
there is no direct evidence for anomalous behavior of the 
interlayer coherence in the BG phase. 

In this Letter, we provide experimental evidence for the 
loss of interlayer coherence at (1/5 ~ l/3)i?$ in the BG 
phase of YBCO by using magneto-optical (MO) imag- 
ing of the critical state field profile. We found that the 
field profile in YBCO with slightly tilted CD's is asym- 
metric, which is explained by the alignment of vortices 
along CD's. The asymmetry, which can be utilized as 
a probe for the interlayer coherence, has a maximum at 
£?$/5 and disappears above B$/3. This result strongly 
suggests that the field-driven boundary exists in YBCO 
in the same field range as BSCCO. 

YBCO single crystals were grown by the flux method 
using gold crucibles |l^]. Rectangular twinned single 
crystals were cut into typical dimensions of 1.0 x 0.5 x 
0.015 mm 3 , so that the edges of the samples are along 
the a and b axes. The critical temperature of the pris- 
tine samples is about 91 K. Crystals were irradiated 
with 600 MeV iodine ions at doses corresponding to 
= 10 kG (crystal A) and B$ = 3 kG (crystal B) 
using the tandem Van de Graaff accelerator with super- 
conducting booster at JAERI. The irradiated direction 
of both samples was tilted 10° from the c axis in the 
y-c plane (see inset of Fig. [l]). Accordingly, we can dis- 
tinguish the effect of CD's from that of twin boundaries 
(TB's). The longitudinal magnetization M parallel to 
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the applied field H was measured by using a commercial 
SQUID magnetometer. We defined Ocd as the angle of 
CD's from the c axis, Oh as the angle of applied field 
from the c axis. The critical state field profile was im- 
aged by using an MO indicator garnet film with in-plane 
magnetization placed on the top surface of the sample in 
the field range ZZ|| C | < 1.5 kOe,\H ±c \ < 1.0 kOe @. 

Inset of Fig. I] shows magnetization hysteresis loop at 
T = 80 K in crystal A with Oh = Ocd- Irreversible 
magnetization shows a maximum at around B§ / 3 and 
we define B p as this peak field. In crystal B, we deter- 
mined B p as a field where [M(0 H = +0 C d) - M(0 H = 
—6cd)]/M(6h = +0cd) shows a maximum, because the 
enhancement of the magnetization by CD's is smaller in 
this sample. Main panel of Fig. [I] demonstrates that in 
both crystals B p /B$ is almost independent of tempera- 
ture and its value is around 1/3 ~ 1/5, which is similar 
to that reported in BSCCO §. 

Figures ||(a) and (b) show typical critical state field 
profiles at H\\ c = 900 Oe in decreasing and increasing 
field branches, respectively. The double-Y shaped cur- 
rent discontinuity lines (ci-lines) are clearly seen, where 
current direction abruptly changes [|I]j2^|| . The cen- 
ter ci-line is significantly shifted in the y direction. The 
shift of the ci-line is much larger than t s tan Ocd , where 
t s is the thickness of the sample. The field dependence of 
the normalized ci-line shift Ay in crystal A is plotted in 
Fig. |^(c). This curve shows hysteresis, which is symmet- 
ric with respect to H = 0. When the field is increased, 
the d-line at the center of the sample shifts towards one 
of the edges (Fig. |2|(b)), whereas it shifts to the opposite 
direction when the field is decreased (Fig. [^(a)). The di- 
rection of the shift is always along the y axis, which is the 
same as that of the inclination of CD's. The hysteresis 
can be summarized as follows. When H and M have the 
same sign, the shift is positive, and if they are opposite, 
Ay is negative. When we reverse the field sweep direc- 
tion, the critical current direction in the sample (or the 
sign of M) is reversed, and then the <i-line shifted to one 
direction disappears and a new ci-line shifted to the op- 
posite direction appears, as shown by the dotted arrows 
in Fig. ||(c) JIU . A small positive shift Ay at H — can 
be explained by the self field trapped in the sample. It 
should be noted that TB's cannot explain the observed 
shift of the ci-line, since TB's run randomly along [110] 
and [110] directions and the symmetry of TB's is different 
from that of the ci-line motion. 

To clarify the relationship between Ay and the align- 
ment of vortices along CD's, we investigated the field 
profile under tilted fields. Figure [3] shows the normalized 
ci-line shift Ay as a function of the misalignment angle 
(Oh — Ocd) in crystal A. A sign change of Ay occurs 
around Oh — Ocd ~ in both increasing and decreasing 
field branches. When \0h — 0cd\ is large, |Ay| becomes 
smaller. This result shows that the misalignment of the 
field from CD's is an important parameter to determine 



the shift of the ci-line. 

Previously, an asymmetric field profile was reported in 
YBCO with tilted CD's from the c axis by Schuster et 
al. pl|] . They observed the in-plane anisotropy of J c , 
and discussed it based on the difference in the nucle- 
ation energy of kinks in two cases, along and across the 
tilted CD's p4| . They considered that the asymmetric 
field profile originates from the different kink nucleation 
between the top and bottom surfaces because of the dif- 
ference of the surface quality. However, we checked that 
Ay at the top and bottom surfaces have opposite polar- 
ities (see Fig. f|(d)), indicating that the asymmetry is a 
rather intrinsic property of vortex systems. In addition, 
we performed MO imaging in the sample with non-tilted 
CD's (Ocd ~ 0°) under tilted fields and confirmed that 
the asymmetric field profile depends only on the misalign- 
ment of CD's and H. 

We propose a model to explain our observations. In our 
model, the shift of the ci-line is caused by the alignment 
of vortices along CD's, which is schematically shown in 
Fig. f|(a). Even when the field is applied away from the 
direction of CD's, vortices can be aligned by CD's if the 
misalignment is not so large. To compensate the differ- 
ence in the directions of H and B, in-plane magnetization 
(My) is induced as shown in Fig. |](a), which is realized 
by the current (J y ) flowing both in and across the Cu02 
planes (Fig. [|(b)). At the same time, out-of-plane mag- 
netization (M z ) is generated by the in-plane current as 
shown in Fig. ^(c). Actual current density in the sam- 
ple is the sum of both currents, and it is limited by J c ±, 
where J c m and J c ± are the in-plane critical current densi- 
ties parallel and perpendicular to the y axis, respectively 
p3j . The existence of J y breaks the balance between J c \ 
and J C 2, because J c \ is always anti-parallel to J C 2- This 
imbalance makes the shift of the ci-line, since the MO 
indicator can detect only the out-of-plane induction B z 
close to the top surface. 

To check this idea, we calculate B z (y) assuming a cur- 
rent density distribution J c ±(y,z) = ±|J c j_ as shown 
in Fig. 4(d), where the contribution of J y is introduced. 
The calculated B z (y) in Fig. 4(e) shows a ci-line whose 
position is shifted from the center of the superconductor, 
just as we observed. One may notice that the in-plane 
magnetization produces small stray fields near the edges 
(y ~ ±1), but in the total B z (y) this effect is negligibly 
small. Our images in Fig. 2 are qualitatively consistent 
with this calculation. 

Our model naturally explains (1) the hysteresis ob- 
served in Fig. ^ by the change of the direction of J c ±, (2) 
the opposite polarity between top and bottom surfaces, 
and (3) the Ay sign change at 0h ~ Ocd by the change of 
the direction of J y . In the full critical state of a rectangu- 
lar sample, we can estimate the relative critical current 
densities in the four regions separated by the double Y- 
shaped ci-lines Q. In Fig. g, the ratio (1 + Ay)/(1 - Ay) 
gives J c i/Jc2, and Ay is therefore equal to the ratio of 
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the current densities J y /J c ±. An important point is that 
the measurements of the shift Ay is much easier than the 
global transverse magnetization (M y ) measurements p3| 
in thin samples, since J y /J c ± can be large due to the 
small thickness even if the integrated M y is small. 

A useful implication of our model is that when vortices 
lose the interlayer coherence with zigzag-like structure 
along z direction, the in-plane magnetization M y and 
hence Ay will disappear. Therefore, the asymmetry of 
the field profile in the critical state can be a powerful 
probe for the interlayer coherence. 

Next let us discuss what happens on the asymmetry 
when we cross B p in sample B. Figure |^ shows J y J J c ± as 
a function of applied field H at T — 55 K and 9h = 0°. 
The lower field part of Fig. | (\H\ < 0.4 kOe) is con- 
sistent with the hysteresis in sample A (Fig. ||(c)). Sur- 
prisingly, the absolute value of J y / J c ± has a maximum 
at ~ Bq> /5 and becomes almost zero (with no hysteresis) 
above Bg,/3. This behavior indicates that vortices are 
aligned along CD's below -B$/3, and the interlayer co- 
herence is suppressed above I?$/3 [|6|. This field range 
of the interlayer coherence anomaly is quite similar to 
that of field-induced recoupling in the VL phase of irra- 
diated BSCCO observed by JPR experiments |J. Fur- 
thermore, at the same field range, the enhancement of J c 
is observed in the BG phase (Fig. |l|), also consistent with 
BSCCO §. This similarity between YBCO and BSCCO 
implies that the field-driven _B$/3 anomaly does not de- 
pend on the anisotropy of the system. Actually, recent 
JPR studies in irradiated BSCCO with different oxygen 
contents have shown no dependence of anomaly field on 
the anisotropy p?| ]. Moreover, our results that the inter- 
layer coherence shows a dramatic decrease at ~ -B$/3, in 
contrast to the increase in the VL state, are consistent 
with the prediction of the simulation result |l3| in the 
BG phase. 

Finally, the reason why this anomaly field is around 
(1/5 ~ 1/3)5$ is still an open question. It is well known 
that the matching effect is observed at H = £?$ |2Sf| 
when the distribution of the pinning center is periodic 
p9| . In the irradiated crystals, however, CD's are ran- 
domly distributed, which suggests that the statistical av- 
eraging may be important to understand the underlying 
mechanism of this number. 

In summary, we observed an asymmetric critical-state 
field profile in YBCO with CD's when the field is tilted 
away from CD's. The asymmetry depends on the field 
sweep direction and the misalignment of the field from 
CD's. We interpret this asymmetry in terms of the in- 
plane magnetization, which is originated from the align- 
ment of vortices along CD's. We proposed that the asym- 
metry of the critical state field profile can be used as a 
powerful probe of the interlayer coherence. The coher- 
ence of vortices along CD's has a maximum at ~ £?$/5 
and becomes small above B$/3. This result in YBCO is 
analogous to the results of JPR in BSCCO. 
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FIG. 1. Characteristic peak field B p normalized by the 
matching field £?$ as a function of temperature (see text). 
(Inset) Magnetization hysteresis loop at T = 80 K in crystal 
A (_B$ = 10 kG). The shaded area shows the field region 
where we nerformed MO imafin?. Arrows indicate the nea.k 
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FIG. 2. Typical critical-state field profiles in the (a) de- 
creasing and (b) increasing field branches at Hu c = 900 Oe 
in crystal A. The bright region corresponds to higher fields. 
The direction of CD's is tilted 10° from the c axis towards 
the y direction. Note the asymmetry of the profile, (c) Hys- 
teresis loop of the normalized d-line shift (Ay) at 80 K in 
crystal A. Inset shows the definition of Ay. The sign of Ay 
is positive when the shift direction is towards y axis. Arrows 
with "inc" and "dec" indicate field sweep directions. When 
the field sweep direction is changed, the old ri-line disappears 
and a new d-line is generated (dotted arrows) |23|. 




FIG. 3. Ay as a function of the misalignment angle 
(9 h — 6c d) for the increasing (circles) and decreasing (tri- 
angles) field branches at H« c — 500 Oe in crystal A. Dashed 
line shows the H direction parallel to the c axis. 
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FIG. 4. (a) In-plane magnetization M y is generated due 
to the misalignment of B and H, when the vortex is par- 
tially trapped. In-plane (b) and out-of-plane (c) components 
of magnetization and their accompanying currents in the field 
decreasing branch at a positive field. Arrows on the sample 
show the direction and magnitude of current density. J c i 
and J C 2 are determined by the constraint that the sum of 
current densities in (b) and (c) cannot exceed the in-plane 
critical current density J c ±. With an assumed current dis- 
tribution (J y /Jc± = 0.29, J _L z,y) in a finite strip sample 
(w : I : d = 675 : 350 : 15), field profile B z (y) along the bro- 
ken line in (d) is calculated (e). Induction profiles from J y , 
Jci+Jc2 and J c x are shown in broken, dotted and solid lines, 
respectively. The divergence at the d-line is clearer in the MO 
image, because the intensity is proportional to B\ . 
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FIG. 5. Jy/J c ± as a function of the applied field H at 
T = 55 K in crystal B (B& = 3 kG). Arrows with "inc" and 
"dec" indicate the field sweep directions. 
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